Accepted for publication in November 10, 2002 edition of the Astrophysical Journal 
Preprint typeset using I^IgX style emulateapj v. 14/09/00 



O 
O 

(N 



O 
m 
o 
o 

(N 
O 

6 



EXTREME ULTRAVIOLET ABSORPTION LINES IN Lya FOREST ABSORBERS AND THE OXYGEN 
ABUNDANCE IN THE INTERGALACTIC MEDIUM 

Randal C. Telfer', Gerard A. Kriss 12 , Wei Zheng 1 , Arthur F. Davidsen 13 , David Tytler 4,5 

Accepted for publication in November 10, 2002 edition of the Astrophysical Journal 

ABSTRACT 

We create stacked composite absorption spectra from Hubble Space Telescope Faint Object Spectrograph data 
from four quasi-stellar objects to search for absorption lines in the extreme ultraviolet wavelength region asso- 
ciated with Lya forest absorbers in the redshift range 1.6 < z < 2.9. We successfully detect O V A630 in Lya 



absorbers throughout the 10 to 10 cm column density range. For a sample of absorbers with 10 < 
A/(HI) < 10 14 2 cm" 2 , corresponding to gas densities ranging from around the universal mean to overdensities 
of a few, we measure an O V A630 equivalent width of 10.9 ± 3.7 mA. We estimate the detection is real with at 
least 99% confidence. We only detect O IV A788, O IV A554, O III A833, and He I A584 in absorbers with Lya 
equivalent widths > 0.6 A, which are likely associated with traditional metal-line systems. We find no evidence 
in any subsamples for absorption from N IV A765, Ne V A568, Ne VI A559, Ne VIII AA770, 780, or Mg X AA610, 



13.2 



< 



625. The measured equivalent widths of O V suggest values of (O V/H I) in the range —1.7 to -0.6 for 10 
A/(H I)< 10 15 cm" 2 . The lack of detectable O IV absorption except in the strongest absorption systems suggests a 
hard ionizing background similar to the standard Haardt & Madau spectrum. Using photoionization models, we 
estimate that the oxygen abundance in the intergalactic medium with respect to the solar value is [O/H] « -2.2 
to —1.3. Comparing to studies of C IV, we estimate [O/C] w 0.3 to 1.2. The overabundance of oxygen relative to 
carbon agrees with other low-metallicity abundance measurements and suggests enrichment of the intergalactic 
medium by Type II supernovae. 

Subject headings: galaxies: intergalactic medium — quasars: absorption lines — ultraviolet: galaxies 



1. INTRODUCTION 

Early studies of the intergalactic medium (IGM) identified 
two separate groups of absorbers (e.g., Sargent et al. 1980): (1) 
"metal-line systems" which have observable metal lines and 
the strongest Lya absorption, typically with a Lya rest-frame 
equivalent width > 1 A, and (2) weaker absorbers having no ob- 
servable metal lines and comprising what is generally thought 
of as the "Lya forest". However, it was suspected that this dis- 
tinction was due at least partly to observational limitations in 
the observability of weak metal lines (Tytler 1987), a suspicion 
that has been verified as technological and computational ad- 
vances have allowed astronomers to detect metal absorption to 
much lower column densities. Constraining the abundances of 
metal ions in the IGM has become an active area of research, 
since this information can in principle be used to determine the 
chemical abundances in the IGM as well as shape and strength 
of the typical metagalactic ionizing background radiation from 
stars and quasi-stellar objects (QSOs). Of particular interest 
is the low column density regime (N(H I)< 10 14 5 cm" 2 ), since 
these systems are presumably far away from local sources of 
metals and ionizing photons and thus offer the best insights into 
enrichment mechanisms and the ionizing background. 

The production of composite or stacked absorption spectra 
of many Lya absorbers is an obvious way to increase obser- 
vational sensitivity. Norris, Peterson, & Hartwick (1983) used 
the composite technique to claim a detection of O VI AA1032, 
1038 absorption. The reality of this detection was challenged 



in a similar study by Williger et al. (1989), but O VI in the IGM 
was detected unambiguously by Lu & Savage (1993) by creat- 
ing a composite absorption spectrum of C IV absorbers. Com- 
posite spectra were used to detect C IV absorption in systems 
withA/(H I)> 10 14 cm" 2 at high redshift (z > 1 .7) by Lu (1991), 
but Tytler & Fan (1994) failed to detect C IV associated with 
weaker lines. In a study similar to Lu (1991) but at z < 0.8, 
Barlow & Tytler (1998) found an order of magnitude stronger 
C IV absorption than Lu (1991). 

A more precise way to increase sensitivity is to build big- 
ger telescopes and better instruments in order to acquire higher 
signal-to-noise ratio (S/N), higher resolution data. With the 
HIRES instrument on the 10m Keck I telescope, individual 
C IV features can be detected down to a column density of 
N(C IV) ~ 10 12 cm" 2 . The majority of Lya absorbers with 
A/(HI)> 3 x 10 14 cm" 2 exhibit such absorption (Tytler et al. 
1995; Cowie et al. 1995; Songaila & Cowie 1996). Besides 
allowing for the direct detection of weaker lines, the existence 
of such high quality data has opened the door to an additional 
detection method introduced by Cowie & Songaila (1998), that 
of comparing the individual pixel optical depths of Lya to 
the expected position of corresponding metal lines. Cowie & 
Songaila (1998) and Ellison et al. (2000) use this technique to 
detect C IV absorption down to an H I optical depth of thi < 1 . 
Schaye et al. (2000) have used the optical depth technique to 
detect O VI to optical depths down to thi ~ 0. 1 . 

Studies have concentrated on C IV not only due to the fact 
that it is expected to be one of the strongest features associated 
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with the IGM (Rauch, Haehnelt, & Steinmetz 1997; Hellsten 
et al. 1998), but also largely because of its near ideal location in 
wavelength space. With a wavelength of 1549 A, this absorp- 
tion feature is near enough to Lya A1216 that both features can 
easily be acquired with the same observations, but far enough 
above Lya that the feature will be longward of the Lya for- 
est for a significant span in redshift. Although C IV is a good 
tracer of metals in the IGM, in standard photoionization models 
the O VI AA1032, 1038 A doublet is expected to be a stronger 
feature at H I column densities < 10 15 cm -2 (Hellsten et al. 
1998). Unfortunately, although the wavelength of C IV makes 
it very appealing for absorption studies, the O VI AA1032, 1038 
A doublet could hardly reside in a worse portion of the spec- 
trum. The region of O VI absorption is muddled not only by 
Lya absorption but by higher order Lyman lines as well, mak- 
ing detections much more difficult. Fortunately, Schaye et al. 
(2000) have overcome this difficulty to detect O VI in the dif- 
fuse IGM, although they make no quantitative claims as to the 
implied amount of O VI in the IGM. 

The extreme ultraviolet (EUV) absorption region has been 
largely ignored, not because this region contains no useful in- 
formation, but because the best available data in this region, that 
from the Hubble Space Telescope (HST) Faint Object Spectro- 
graph (FOS), is inferior in resolution and S/N to that used at 
longer wavelengths. However, some of the strongest expected 
absorption features occur in this region, in particular O V A630 
and O IV A788 (Verner, Tytler, & Barthel 1994). The O V fea- 
ture is especially worthy of interest because its large oscillator 
strength (/ = 0.514, greater than H I Lya) should make it com- 
petitive in strength to the O VI feature even at column densities 
down to ~ 10 13 cm -2 . Despite the fact that the FOS has rather 
low resolution by modern standards (R m 1300), searching for 
O V rather than O VI has the advantage that the features are in 
a spectral region where the Lya forest is less dense owing to 
the strong evolution of the density of lines with redshift. The 
detection of oxygen lines in the IGM can be used to infer the 
oxygen abundance, which is expected to be overabundant rel- 
ative to carbon with respect to the solar ratio if the early uni- 
verse is enriched by Type II supernovae. Strong evidence for 
such abundance patterns are seen at high redshift in Lyman- 
limit systems (Reimers et al. 1992; Kohleret al. 1999), damped 
Lya absorbers (Pettini, Lipman, & Hunstead 1995; Lu et al. 
1996), and C IV absorbers (Dave et al. 1998; Songaila 1998). 
Here we extend this work to the lower column density regime 
by carrying out a statistical search for absorption in the EUV 
using the presently available data. 

In §2 we discuss the data and the spectral selection crite- 
ria that we use for this study. We begin §3 by discussing the 
characterization of our Lya absorber sample from optical data. 
This is followed by a presentation of the details of the technique 
that we use to generate the stacked absorption spectra and the 
method for measuring the equivalent width of the resulting ab- 
sorption features. In §4 we describe the results of our search 
for various samples selected by the strength of Lya absorption, 
particularly the detections of O V. We consider the implications 
of our results in §5 by comparing the measurements and limits 
on the equivalent widths of the absorption features to simulated 
data to constrain the abundance ratios of the relevant ions to H I. 
Using photoionization models, we then infer the abundance of 
oxygen in IGM. We provide a brief summary of our results in 
§6. 



2. DATA 

We have searched the HST archives for QSO spectra suitable 
for detecting EUV absorption features in the Lya forest, with 
the primary goal being the detection of O V A630. To include 
an object in this study, we require that: 

1 . Data are available shortward of 630 A in the QSO 
emission frame, making the O V feature accessible. 

2. The data are obtained with the FOS high-resolution 
gratings. 

3. The S/N in the region of the O V absorption is high 
enough that the rms deviation of the pixels is dominated 
by absorption features rather than Poisson noise. 

For the typical density of lines in the applicable 
redshift range (1.6 < z < 2.9), the necessary S/N is 
~ 5 per pixel. At this S/N, we can make a fairly 
reliable determination of the continuum level, which is 
necessary for our study. 

4. We have access to high-quality optical spectra for 
determining the redshift and strength of the Lya forest 
lines. 

We have found four QSOs that match these criteria. The FOS 
data for these QSOs were reduced in 2001 July with the best 
available reference files. In Table 1 we list these objects, 
their redshifts, the total exposure time of the FOS data with 
each grating, and the dates the data were obtained. The zero- 
point shift for the wavelength scale given in of the last column 
will be discussed in §3.3. All of the data were obtained with 
the red Digicon detector on the FOS. A G190H spectrum of 
HS 1700+6416 was obtained with the blue Digicon on 1997 
May 29. However, the data are of somewhat lower resolution 
and for consistency we use only the data from the red Digicon. 

For all four objects we have access to Keck HIRES spectra 
for examining the H I Lya forest. The spectrum for HE 2347- 
4342 is courtesy of A. Songaila and is discussed in Songaila 
(1998). Two of the other three spectra are discussed elsewhere, 
and we point the reader to the corresponding articles for details: 
Kohler et al. (1999) for HS 1 103+6416 (1997 April 9 data only) 
and de la Varga et al. (2000) for HE 1 122-1649. The spectrum 
of HS 1700+6416 consists of a 4300 second exposure on 1995 
May 10. This spectrum covers wavelengths from 3727 A to 
5523 A with a typical S/N of ~ 50 per resolution element. 

3. ANALYSIS 

3.1. Lya Forest Data 

The Keck HIRES data are normalized to the continuum by 
fitting a spline by hand, placing nodes at what appear to be un- 
absorbed portions of the spectra. We then fit the spectra in sec- 
tions of <~ 50 A down to the redshifted position of Ly/3. We use 
the IRAF task specfit to fit the individual absorption features 
with Gaussians in optical depth: 



r = ro exp 



A-Aq 
X b/c 



(1) 



Each absorption feature is thus described by three free parame- 
ters: the central wavelength Ao, the optical depth at line center 
To, and the Doppler width b. In a few cases where the blending 
is severe, the ^-parameter is fixed at 30 km s" 1 , a typical value 
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Table 1 
FOS Data 



Object 


z 


Grating 


Exposure (s) 


Date 


Shift (A) 


HS 1103+6416 


2.190 


G270H 


5336 


1996 Oct 31 


+0.10 






G190H 


8628 


1996 Oct 31 


+0.59 


HE 1122-1649 


2.400 


G270H 


780 


1996 Jul 18 


+0.30 






G190H 


5620 


1996 Jul 18 


+1.25 


HS 1700+6416 


2.743 


G270H 


9000 


1992 Dec 13 


0.00 






G190H 


16120 


1992 Dec 12-13 


+0.85 


HE 2347-4342 


2.885 


G270H 


1320 


1996 Jun 7 


+0.80 






G190H 


7930 


1996 Jun 7 


+1.05 



for Lya forest lines (Rauch et al. 1992; Hu et al. 1995), to pre- 
vent b from taking on an unreasonable value. We fit all lines 
with peak optical depths > 0.1, although lines weaker than this 
were sometimes fit to deblend stronger lines. 

We impose several selection criteria on our list of Lya forest 
lines: 

1 . As is frequently done in Lya forest studies, we ignore 
lines shortward of Ly/3 in the rest frame of the QSO to 
avoid confusion of Lya lines with higher-order lines. 

2. We discard lines within 5000 km s" 1 of the QSO 
rest frame in order to avoid contamination from the 
proximity effect (see e.g. Scott et al. 2000). This should 
also eliminate many associated absorbers, which often 
occur near the QSO in velocity space. However, some 
known associated absorbers are displaced substantially 
more than 5000 km s" 1 blueward of the QSO redshift 
(Jannuzi et al. 1996; Hamann, Barlow, & Junkkarinen 
1997). We assume that such systems are rare and not a 
significant source of contamination in our sample. 

3. We do not use lines that correspond to known 
Lyman-limit systems. These high-column-density 
systems have been well-studied on an individual basis 
(Vogel & Reimers 1995). Since we have data at the 
Lyman limit for all candidate lines by virtue of our 
choice of sample, it is easy to identify Lyman limit 
features with r > 0.1, corresponding to H I column 
densities > 10 16,2 cm" 2 . 

4. We eliminate lines that are identified as likely metal 
lines associated with Lyman-limit or metal-line 
absorption systems. 

5. We ignore features in wavelength regions that are 
contaminated by strong clusters of such metal lines. 

6. We dismiss lines with b < 14 km s" 1 as probable 
unidentified metal lines. Detailed studies of the 
distribution of b (Hu et al. 1995; Kirkman & Tytler 
1997) show that there is a sharp cutoff at around 20 
km s" 1 and virtually no lines with b < 14 km s" 1 . 

7. We also exclude lines with b > 80 km s" 1 . Although a 
few of these lines are probably real, they are certainly 
atypical. Some are likely due to blending. 

It is straightforward to use the best-fit Gaussian parameters 
to calculate the equivalent width, redshift, and column density 



of H I. To illustrate that we have extracted a reasonable distri- 
bution of lines, we plot the number of lines in our sample as a 
function of column density in Figure 1 . The median redshift of 
the lines in our sample is z = 2.35. We have placed all of the 
lines between 10 14 2 and 10 16 2 cm" 2 into a single bin since the 
column densities of these strongly saturated lines are poorly de- 
termined. In principle, we could use higher-order Lyman lines 
to measure the column density. Unfortunately, for most of the 
Lya features in our sample, the corresponding Ly/3 feature is 
either in the low wavelength portion of the Keck data, where 
the S/N is low, or at wavelengths in between the Keck and FOS 
spectral coverage where we have no data. 
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FIG . 1 . — Distribution of the number of lines in our sample with column den- 
sity. The dashed line is our best fit to the distribution between 10 13 and 10 14,2 
cnT 2 with /3 = 1.41 ±0.05. The dotted line represents [3 = 1.83 connecting to 
our best fit at 10 14 ' 3 cm 4 . 

The ordinate in Figure 1 is proportional to the column den- 
sity distribution function, often called f(N) in the literature. As 
is standard practice, we fit the distribution between 10 130 and 
10 14 - 2 cm" 2 with a power law, f(N) oc N' 13 . Our best-fit power- 
law index is (3 = 1.41 ±0.05, plotted as the dotted line in Fig- 
ure 1, in agreement with Kim et al. (1997) for this range of 
column densities at z « 2.3. The fit appears to be good down 
to 10 13 cm" 2 , which argues that we are essentially complete to 
this column density. The column density distribution function 
is known to steepen above 10 143 cm" 2 (Kim et al. 1997). We 
adopt [5 = 1.83 as found by Petitjean et al. (1993) to be a good 
fit between - 10 14 and 10 16 cm . We connect this distribution 
to our best fit at 10 14,3 cm" 2 and plot it as the dotted line in Fig- 
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ure 1 . Clearly this agrees well with our high-column-density 
point. When considering the distribution of column densities 
for 7V(H I)> 10 14 3 cm" 2 which is important for our data simu- 
lations in §5.3, we therefore adopt [3 = 1 .83 as representative of 
our data. 



3.2. Combination Technique 

The process of creating composite absorption spectra is con- 
ceptually straightforward. First we combine the G190H and 
G270H spectra into a single spectrum so that the wavelengths 
in the overlap region do not receive double weighting. To facil- 
itate this, we resample the individual spectra into 0.1 A bins in 
the observer frame - a factor of several smaller than the original 
pixels so as not to alter the character of the data. Then we nor- 
malize the FOS spectra to the continuum using a cubic spline 
in the same way as we did for the optical data. Instrumental 
artifacts and a damped Lya absorption feature in one object 
are masked out. We also remove the noisy portion at the short 
wavelength end of the G190H data, below ~ 1650 A, where the 
sensitivity of the red Digicon detector becomes low and the S/N 
drops per original FOS pixel drops below 5. 

For each Lya forest feature in the particular sample, the com- 
bined FOS spectrum of that object is shifted into the rest frame 
of the absorber by dividing the wavelength scale by 1 +z a bs- 
These shifted spectra are then rebinned into a common wave- 
length scale in the absorber frame for coadding. For conve- 
nience we choose a bin size of 0.1 A, which corresponds to 
48 km s" 1 in the region of O V A630, slightly smaller than the 
typical pixels in the original data. 

For each pixel in the proposed composite spectrum, we then 
have a distribution of pixel fluxes from the individual shifted 
spectra. A crucial question remains - what statistic of this 
distribution do we use as the composite flux in order to opti- 
mize the sensitivity of our measurements? The answer depends 
strongly on the character of the noise, which in our study is 
due largely to random absorption features. If we were dealing 
with high S/N data and the continuum in the spectral region 
of interest were not polluted by random absorption lines, as is 
typical in searches for weak C IV absorption, then the noise 
would be Gaussian distributed and a mean or median combina- 
tion method would do nicely to minimize the noise in the result. 
However, when the features of interest are strongly blended into 
the Lya forest the noise is distinctly non-Gaussian. 

On the top of Figure 2 we plot as a solid line the histogram of 
the pixel fluxes of the FOS spectra of our four QSOs, after the 
G190H and G270H data have been binned to 0.1 A and com- 
bined. Clearly the pixels have a very asymmetric distribution 
in flux with a long tail to low flux due to the many absorption 
lines. We expect that at any random wavelength, where there is 
no feature, the contribution to the composite from the individ- 
ual shifted spectra will consist of a random sampling from this 
distribution. We would like to maximize the S/N in these ran- 
dom samplings, which for our composite spectra corresponds to 
the S/N in the continuum. We can take any general percentile 
of the data by linearly interpolating the cumulative distribution, 
where the lowest flux at a given pixel is defined as the zeroth 
percentile and the highest flux as the 100th. The simplest and 
most common example is the median or 50th percentile, but we 
will show that we can increase the S/N significantly by making 
a different choice. 
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FIG. 2. — Top: Distribution of individual pixel fluxes in the FOS data of the 
four QSOs, after the spectra have been normalized and resampled into 0.1 A 
bins. Bottom: The same distribution displayed as a cumulative distribution. 

The optimal percentile should be near the peak of the dis- 
tribution, which one can see from the cumulative distribution 
in Figure 2 is around the 80th percentile. We perform Monte 
Carlo simulations to show explicitly that this is the case. We 
take 100000 random samplings of 100 points from the flux dis- 
tribution, each time calculating the nth percentile of the sam- 
ple. Figure 3 shows the result of the simulations as a function 
of percentile, where the S/N is defined as the average value of 
the nth percentile of the sample divided by the rms deviation. 
Indeed, we see that the S/N peaks where the flux distribution 
peaks, with the optimal choice being the 79th percentile. An 
earlier version of these simulations led us to conclude that the 
78th percentile was optimal which was used in the subsequent 
analysis. Since Figure 3 shows that the difference is negligible 
(a 0.3% difference in S/N), we use 78th percentile spectra for 
this work. The 78th percentile yields a S/N nearly 40% greater 
than the median in these simulations. The mean yields a S/N of 
49. 1, similar to the median. The optimal percentile is not a gen- 
eral result, but rather is peculiar to the S/N and absorption line 
density of our data. This discussion gives our method a solid 
statistical basis and points out clearly the importance of care- 
fully considering the noise distribution of a particular dataset in 
order to maximize the resulting S/N in the composite. In §4 we 
will demonstrate empirically that taking the 78th percentile of 
the distribution is indeed essentially optimal for our data. 

Since we have selected our QSO sample to search for O V 
A630, we have FOS data at this wavelength for virtually all of 
the Lya absorbers on our sample, the only exceptions being 
those few for which the wavelength of the O V absorption hap- 
pens to lie in a masked region. Overall, we have FOS data for 
over 90% of the absorbers for rest-frame wavelengths between 
560 and 890 A. At wavelengths longer that ~ 900 A, the S/N 
in the composites becomes very low due to the high density of 
Lyman lines. As a result, we cannot place any interesting con- 
straints on metal lines at these wavelengths. There are some 
interesting lines shortward of 560 A, including O IV AA553, 
554 and Ne VI A559. Although we could use this O IV doublet 
to constrain the amount of O IV, we concentrate on O IV A788 
because it provides comparable sensitivity and we have data for 
more of the absorbers. For all other ions that we consider there 
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is a single obvious choice for which line to use. 




40 60 
Percentile 

FIG. 3. — Average S/N of 100000 random samplings of 100 points from the 
distribution of all the pixel fluxes from the FOS spectra shown in Figure 2, as 
a function of the percentile taken for each sample. 



3.3. Wavelength Calibration 

There is a well-documented problem with the absolute wave- 
length scale of the FOS (Rosa et al. 1998). A solution to this 
problem has been found for data from the blue detector, which 
has a systematic shift over time, but a fix for the red detector 
data has not yet been established. We determine this empiri- 
cally for our data. First, we find the relative shift between the 
G190H and G270H data for each QSO. Using data in the region 
of wavelength overlap from around 2220-23 10 A, we calculate, 
to the nearest 0.01 A, the wavelength shift which minimizes the 
rms deviation of the pixels. We then shift the G190H data by 
the required amount and combine the two into a single spec- 
trum. This method works well and gives shifts consistent with 
what we estimate by measuring individual features in the over- 
lap region. 

The overall shift of the combined G190H / G270H data on 
an absolute scale is more difficult to determine. Given the S/N 
of the data, Galactic absorption lines are not readily identifi- 
able for HE 1122-1649 and HE 2347^1342, and only a few 
for HS 1103+6416 and HS 1700+6416, so in general using 
the Galactic lines is not a reliable way to set the wavelength 
scale. To estimate the overall shift, we create composite ab- 
sorption spectra of the individual objects using only Lyct forest 
lines with Lyct equivalent widths > 0.5 A. The lower limit was 
chosen so that a well-defined O V line could be seen for each 
object. In each case, a strong absorption feature occurs near 
the expected position of O V A630, but slightly shifted by up 
to ~ 100 kms" 1 . We shift the FOS spectra by the amounts 
listed in Table 1 to center the O V feature as near as possible to 
629.73 A. The shifts listed for the G190H data include both the 
overall shift and that necessary for agreement with the G270H 
data. The required shifts are modest and well within the errors 
of the FOS red detector. We use these shifted spectra for all of 
our final results. To check the accuracy of the applied shifts, 
we generate a composite using absorbers with Lyct equivalent 
widths > 0.6 A for each object. The larger lower limit makes 
some lower-ionization lines more visible, as discussed in §4.3. 
For each object, we have measured the positions of other lines 



for which we have data and that are strong enough to make an 
accurate determination of the wavelength centroids. For this 
purpose we use lines from the following list: Ly/3, Ly7, Ly<5, 
He I A584, O IV A554, and O IV A788. We find that the aver- 
age shift of the lines for each object is < 40 km s" 1 . Thus, the 
residual uncertainty in the wavelength scale is significantly less 
than the instrumental resolution and should have little effect on 
our results. 

One may argue that by applying these shifts we could be ar- 
tificially enhancing the O V feature, since we have used O V 
to set the wavelength scale. This may be so, but we point out 
that such an argument would apply only to O V and only for 
absorbers with a Lya equivalent width > 0.5 A. Our important 
results are derived almost entirely from weaker lines, and there 
is very little overlap between the relatively weak absorbers we 
analyze in detail and the strong systems we use to set the wave- 
length scale. 

3.4. Measuring the Equivalent Width 

To measure the equivalent widths of the composite absorp- 
tion features we use the profile weighting method as described 
by Lu (1991). Since this is a crucial part of our analysis we 
describe it briefly here. For an unresolved absorption feature, 
if we know the instrumental line profile and the central wave- 
length of the feature, then in principle every pixel gives an 
estimate of the absorption equivalent width. To get an accu- 
rate measurement of the equivalent width, we should take a 
weighted mean of the individual pixel estimates of the equiv- 
alent width. If the error is constant as a function of wavelength, 
which is roughly true for our data, the optimal weighting is the 
square of the instrumental profile. The equivalent width is then 
calculated as 



EW = 



(2) 



where /} is the flux, c, is the continuum level, and is the 
normalized instrumental profile. For the instrumental profile 
we use a Gaussian with a full-width at half maximum derived 
from the composite Ly/3 features, which we find to be typically 
around 290 km s" 1 . For the continuum determination we use a 
simple linear regression fit to the flux between ±2500 km s" 1 , 
eliminating the region between ±500 km s" 1 where the feature 
lies. 

This weighting technique yields accurate results by relying 
on the fact that every individual metal feature associated with 
the sample has not only the same observed profile shape, but 
also, for any given ion, the same wavelength centroid in the ref- 
erence frame of the absorber. For our sample, we determine 
the redshift of the absorption systems from the Lya absorption 
lines. As discussed in §3.3, the mean shift of absorption lines 
with respect to this reference frame due to uncertainties in the 
FOS wavelength scale is < 40 km s" 1 for all four lines of sight. 
However, what are the uncertainties for individual absorbers? 
Ellison et al. (1999) find an rms dispersion in the velocity off- 
set between C IV and H I absorption of 27 km s" 1 , while Ellison 
et al. (2000) find a dispersion of 17 km s" 1 . The authors point 
out that these offsets can have a dramatic effect on the results 
of the stacking technique when applied to high-resolution data, 
such as that from the Keck HIRES. However, at the FOS reso- 
lution, typical offsets of 17-27 km s" 1 should have little effect 
on the composites. 
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FIG. 4. — Composite spectra of Lya forest features with N(U I)< 10 14 2 cm 2 for various minimum H I column densities, plotted in the velocity space of O V 
A630. 



4. RESULTS 

We divide our analysis into three groups based on the Lya 
feature strength. One group consists of the Lya lines with 
N(H I)< 10 14 2 cm" 2 . These lines are unsaturated and therefore 
we have accurately measured column densities, so the analy- 
sis requires fewer assumptions than for saturated lines. This 
is the most important group because it samples absorbers for 
which little information is available about metal content. For 
stronger lines, we cannot reliably measure the column density 
of H I, so we must instead discriminate absorbers by Lya equiv- 
alent width. We divide the stronger lines into two groups - an 
intermediate sample with Lya equivalent widths < 0.6 A, and 
a sample with Lya equivalent widths > 0.6 A. Although the 
exact division is somewhat arbitrary, this strongest group con- 
sists of lines with H I column densities > 10 15 cm" 2 , which are 
generally associated with metal-line systems. Each group thus 
spans roughly an order of magnitude in column density. We will 
discuss the low-column-density sample first, since these weak 
lines are useful for assessing the success of our composite tech- 
nique. 

4.1. Low Column Densities: N(H /)< 10 14 ' 2 cm" 2 

It is not obvious a priori at what H I column density we 
should place the lower limit of our sample to achieve the most 
secure detection, and in fact this depends on the unknown 
strength of the feature. One would like to include many lines 
to increase the S/N, but this involves adding in lines of lower 
and lower column density. At some point, the feature will be- 
come so weak that including more lines essentially adds only 
noise, thereby reducing the significance of any detection. As 
pointed out by Lu (1991), there should therefore be some opti- 



mal lower limit which maximizes the S/N of the result. To find 
the optimal detection, we let the lower limit range from 10 13 
to 10 13 9 cm" 2 . The resulting 78th percentile spectra are shown 
in Figure 4 in the velocity space of the O V A630 feature. The 
number of lines and the O V equivalent width for each sample, 
measured as described in §3.4, are listed in Table 2. One can 
clearly see how the noise is reduced by including more lines, 
and the O V feature becomes well-defined. In none of the sam- 
ples do we find any evidence for O IV A788. 

In order to assess the errors on these results and their sig- 
nificance, we perform Monte Carlo simulations in which we 
randomize the redshifts of the Lya features so that they can 
take on any value in the allowable range, i.e., between Ly/3 and 
5000 km s" 1 blueward of Lya for the particular object, as de- 
fined previously. The resulting composites should thus consist 
of random noise, and by measuring the equivalent widths in the 
simulated spectra in the same way as for the real data, we can 
evaluate the errors on our equivalent width measurements. We 
have performed 1000 simulations for each sample, and the rms 
deviations of the equivalent width in these Monte Carlo spectra 
are shown as the errors on the O V equivalent width in Table 2. 
The mean equivalent width in these simulations is slightly less 
than zero, typically by -0.5 mA, indicating that there is a slight 
zero-point offset associated with our measurement technique 
that is likely due to the continuum placement. Because for our 
analysis we compare our results to simulated data, which also 
include this effect, we do not add it into our measurement of the 
line equivalent width. 

These simulations can also be used to answer the question of 
how likely it is that the feature is due to random chance. In col- 
umn 4 of Table 2, we list the fraction of simulated equivalent 
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widths that were less than that measured in the real data, which 
can be interpreted as the probability that the result is not due to 
chance. For a lower limit of 10 13 5 cm" 2 , none of the simula- 
tions exceeded the measured value, so we list the probability as 
> 0.999. The significance maximizes with a lower limit around 
j q13.2 j q jo 13 7 cm" 2 . Because so few (two or fewer) of the sim- 
ulations exceed the real measurements for lower limits of 10 132 , 
10 13 - 3 , and 10 13 5 cm" 2 , we cannot reliably distinguish which is 
of greater significance. However, we find that the distribution 
of equivalent widths in the simulations are well-described by a 
Gaussian, and therefore one can simply use Gaussian statistics 
to estimate the significance. The best detections correspond to 
approximately a 3a fluctuation above the mean noise level and 
are significant at greater than a 99% confidence level. 

Table 2 

O V EWs for Weak Lyo absorbers: N(H I)< 10 14 2 crrf 2 



log/V min (Hl) n OVEW(mA) P rea i 



13.9 


69 


19.0±8.9 


0.976 


13.8 


96 


15.7±7.4 


0.987 


13.7 


136 


17.3±6.2 


0.995 


13.6 


163 


14.2±5.5 


0.996 


13.5 


203 


14.9±5.0 


> 0.999 


13.4 


243 


10.8±4.8 


0.987 


13.3 


298 


12.5±4.1 


0.998 


13.2 


347 


10.9±3.7 


0.999 


13.1 


394 


8.5 ±3.7 


0.987 


13.0 


444 


7.3±3.3 


0.956 



Table 3 

Effect of Percentile on O V Feature Significance 



Percentile 


O V EW (mA) 


^real 





132.4 ±148.9 


0.778 


10 


-0.4 ±24.7 


0.568 


20 


7.0±13.3 


0.820 


30 


12.5 ±11.3 


0.904 


40 


13.0±8.7 


0.950 


50 


11.1 ±7.2 


0.957 


60 


10.6±5.9 


0.970 


70 


11.7±5.0 


0.995 


80 


9.7±3.6 


0.997 


90 


7.9±3.8 


0.986 


100 


14.6 ±16.5 


0.813 



Since the low-column-density samples contain many lines, 
we can use them to illustrate the validity of our choice of taking 
the 78th percentile in forming the composite spectrum. As our 
test case we use the 7V(H I)> 10 13 2 cm" 2 sample. At the top of 
Figure 5, we show the 78th percentile spectrum with 347 lines 
in the sample. Below this, we show various other percentiles 
of the contributing spectra, from the 10th (lowest) to the 90th 
(highest), plotted at true scale. We can measure the equiva- 
lent width of the feature in each percentile and calculate the 
noise with Monte Carlo simulations as before. These results are 
shown in Table 3. As expected, although the equivalent width 
of the feature is roughly constant, the noise varies considerably. 
The gain in S/N by using the 78th percentile over, for example, 
the median is even greater than what would be expected based 
on our simulations in §3.2. This is probably because the noise 
is due mostly to random absorption features which span several 



pixels and thus cause correlated errors in the continuum. We 
list the significance of the detection in each percentile in Ta- 
ble 3, and none of these exceed the significance achieved with 
the 78th percentile, although the results are similar for the 80th 
percentile. Based on the simulations in §3.2 and the success of 
the 78th percentile when applied to the real data, we feel justi- 
fied in using this technique. 




I I I I I I I L 

-3000 3000 

Velocity (km s -1 ) 



FIG. 5. — Top: Composite 78th percentile spectrum of all lines with 
10 13 - 2 <N(H I)< 10 14 2 cm" 2 plotted in the velocity space of O V A630, as 
shown in Figure 4. Bottom: Various percentile spectra of the same data show- 
ing the effect of the percentile on the result. The spectra are in steps of 10 
percent, ranging from the 10th percentile (bottom) to the 90th percentile (top). 



4.2. Intermediate Column Densities: 10 14 ' 2 <N(H /)< 10 15 

cm z 

The division between the Lya forest and metal-line systems 
is not a precise concept. In examining the intermediate-column- 
density lines we therefore try a range of upper limits on the Lya 
equivalent widths. The measured equivalent widths of O V are 
listed in Table 4. The spectra are shown in Fig 6, where we 
display not only O V but also O IV A788 to illustrate that there 
is no evidence for a detection, which we quantify in §4.4. As 
opposed to the way we analyzed the low-column-density lines, 
here as we include more lines we are including lines of increas- 
ing strength and thus we definitely expect the significance of 
the detection to increase. This is clearly the case, and with a 
Lya EW upper limit of 0.6 A we can detect O V with > 4er 
confidence. 
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4.3. High Column Densities: N(H I)> 10 15 cm" 2 

For a ^-parameter of 30 km s" 1 , a Lya EW of 0.6 A corre- 
sponds to an H I column density of 10 16 cm" 2 . Since we ex- 
cluded systems with column densities > 10 16,2 cm" 2 based on 
Lyman-limit detections, lines with a Lya EW > 0.6 A must be 
high-column-density systems with large ^-parameters, and we 
expect them to be associated with traditional "metal-line sys- 
tems". Although these strong systems are scarce and not the 
primary interest of this work, we can detect many absorption 
lines at such large column densities which further illustrates the 
success of our technique. 

In Figures 7-9 we show plots for O V A630, O IV A788, O IV 
A554, O IV A608, O III A833, and He I A584. These lines have 
all been observed on an individual basis in Lyman-limit systems 
(Reimers & Vogel 1993; Vogel & Reimers 1993) and metal- 
line systems (Vogel & Reimers 1995). The measured equiva- 
lent widths with errors from Monte Carlo simulations are listed 
in Table 5. The O IV line at 554.07 A has a weaker doublet 
companion at 553.33 A which is visible in Figure 8. We mea- 
sure only the equivalent width of the A554 line, although the 
presence of the A553 line probably affects this measurement 
slightly. All of the lines are detected with fairly high signifi- 
cance (> 0.98) in the EW> 0.6-0.8 A samples, with the excep- 
tion of O IV A608. By itself the presence of this line would be 
tenuous at best, but given the strength of the other O IV lines it 
is present at about the expected level. 



Table 4 

O V EWs for Intermediate Lyq Absorbers: 10 142 < 
N(H I)< 10 15 cnT 2 



Max Lya EW (A) 


n 


O V EW (mA) 


Preal 


0.4 


20 


37.8± 19.1 


0.962 


0.5 


36 


48.3 ±13.9 


0.998 


0.6 


49 


45.3 ±11.4 


> 0.999 


0.7 


55 


48.2 ±10.5 


> 0.999 



4.4. Limits on Other Ions 

Using the same Monte Carlo simulations that we used to de- 
termine error bars on the observed features, we can define 90% 
confidence limits on the equivalent widths of absorption for un- 
observed ions. We do this not only for O IV A788 absorption 
but also for several other ions which have transitions which one 
might expect to observe in the EUV. Specifically, we constrain 
He I A584, O III A833, N IV A765, Ne V A568, Ne VI A559, 
Ne VIII A770, and Mg X A610. We concentrate on three spe- 
cific samples: (1) the absorbers with 10 13 ' 2 <V(HI)< 10 14 ' 2 
cm" 2 , which we refer to as sample L or the low-column-density 
sample, (2) the absorbers with N(H I)> 10 14 2 cm" 2 and Lya 
equivalent widths < 0.6 A, which we will to as sample I or 
the intermediate-column-density sample, and (3) the absorbers 
with Lya equivalent widths > 0.6 A, which we refer to as 
sample H or the high-column-density sample. The equivalent 
width limits for various unobserved ions are listed in Table 6 for 
all three samples. For completeness we include the measured 
equivalent width for detected ions. The ratios of the these ions 
to H I are discussed in §5. 



5. DISCUSSION 

5.1. The O V/HlRatio 

The interpretation of composite spectra is not straightfor- 
ward. In order to reliably interpret the measured equivalent 
widths, it is necessary to generate simulated data to compare 
to the results. As a simplification, we first will assume that the 
mean ratio (O V/H I) is constant within each of our subsamples. 
This may or may not be a good assumption, depending on the 
ionization model, although our subsamples span a small enough 
range of N(H I) that we do not expect (O V/H I) to vary by more 
that a factor of ~ 2. Making this assumption allows us to derive 
results that are, as much as possible, model independent. Later 
in this section we will drop this restriction and simulate data us- 
ing particular ionization models. Our Monte Carlo simulations 
proceed as follows: 

1 . We randomize the redshift of each line to take on any 
value in the allowable range, as we did for our noise 
simulations in §4. 

2. We calculate the O V column density for each absorber 
from the H I column density and the given relative 
abundance (O V/H I). 

3. The equivalent width of O V A630 is calculated from the 
O V column density assuming a particular /?-parameter. 

4. For the particular (randomized) redshift of each 
absorber, we insert a Gaussian absorption feature with 
a FWHM of 290 km s" 1 and the appropriate equivalent 
width at the position of the O V A630 absorption. 

5. The absorbers are combined into a composite and the 
equivalent width of O V A630 is measured as we have 
done for the real data. 

An important assumption in this process is the ^-parameter 
used to calculate the equivalent width. For the weak lines in 
the low-column-density sample, this is fairly unimportant, since 
most of the lines will be reasonably near the linear portion of the 
curve of growth. It will have some effect, though, and it is worth 
considering since this will have a more pronounced effect when 
we model lines at higher column densities. Individual metal 
lines as observed in C IV absorption typically have quite small 
^-parameters. As examples, Cowie et al. (1995) find a median 
b for C IV lines of 10 km s" 1 , while Ellison et al. (2000) find a 
median b of 13 km s" 1 . However, when the C IV absorption is 
strong, the absorption corresponding to particular Lya features 
often consists not of a single line but several such narrow ab- 
sorption features. The effective ^-parameter for such metal-line 
clusters, when viewed as a single absorption feature, is signif- 
icantly larger than 10 km s" 1 . Because individual metal lines 
are not generally detected at the small H I column densities that 
we probe, we do not know the typical structure of the metal- 
line absorption. Strong C IV absorbers often contain perhaps 
3-5 components (see, for example, the C IV profiles in Ellison 
et al. (1999)). An effective ^-parameter of 50 km s" 1 should 
therefore provide a reasonable extreme value for the typical b- 
parameter of the metal absorption associated with a sample of 
H I lines. To illustrate the effect of b on our results, we create 
simulated data for ^-parameters of 10 and 50 km s" 1 . A physi- 
cally realistic picture is probably somewhere in between, since 
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FIG. 6. — Composite spectra of Lya forest features with N(U I)> 10 14 2 cm 2 for several maximum Lya equivalent widths. The region for O V A630 is plotted 
on the left, while the corresponding O IV A788 region is shown on the right. The O IV A788 feature is not detected in any of these spectra. 
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FIG. 7. — Composite spectra of the strongest Lya forest features in our sample for several minimum Lya equivalent widths. O V A630 is shown on the left, O IV 
A788 on the right. 
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FIG. 8. — Same composites as Figure 7, except O IV A554 is shown on the left and O IV A608 on the right. The O IV spectra are plotted in the velocity space of 
the 554.07 A line. Above the O IV spectra we indicate the positions of this line and the doublet companion at 553.33 A. 
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Table 5 

O V, O IV, O III, and He I EWs for Strong Lyo Absorbers: N(H I)> 10 15 cm" 2 



Min Lya EW (A) 








Equivalent Width (mA) 






n 


O V A630 


O IV A788 


O IV A554 


O IV A608 


O III A833 


He I A584 


0.8 


11 


156.6 ±25.6 


97.3 ±39.2 


124.3 ±34.6 


52.3 ±28.6 


120.9±45.2 


110.3 ±29.7 


0.7 


14 


146.1 ±21.7 


91.2±36.7 


68.6±28.7 


30.4 ±25.2 


113.8±40.4 


111.0±25.1 


0.6 


20 


99.6±18.4 


92.2 ±27.7 


77.3 ±22.7 


32.9±21.3 


62.7 ±30.5 


64.0 ±20.5 


0.5 


33 


85.7±13.7 


44.8 ±21.5 


49.2 ±17.2 


19.0± 15.6 


29.6±22.9 


4.7 ±15.4 



Table 6 

Equivalent Widths and Column Density Ratios to H I 





A (A) 


Sample L 


Sample I 


Sample H 


Ion 


EW(mA) 


(Xi/Hl) 


EW(mA) 


(Xi/Hl) 


EW(mA) 


He I 


584.33 


<6.0 


<-1.0 


< 16.4 


<-1.5 


64.0 ±20.5 


NIV 


765.14 


<6.2 


<-1.5 


<23.3 


<-1.9 


< 37.0 


OIII 


832.93 


<8.1 


<-1.0 


<23.7 


<-1.4 


62.7 ±30.5 


oiv 


787.71 


<7.0 


<-1.0 


< 18.8 


<-1.5 


92.2 ±27.7 


ov 


629.73 


10.9 ±3.7 


-1.6 to -0.6 


45.3 ±11.4 


-1.7 to -0.8 


99.6±18.4 


Ne V 


568.42 


<5.1 


<-0.7 


< 16.6 


<-1.0 


<26.3 


Ne VI 


558.59 


<7.9 


<-0.7 


<20.1 


<-l.l 


< 31.2 


Ne VIII 


770.40 


<8.2 


<-0.8 


<25.3 


<-l.l 


< 38.9 


MgX 


609.79 


<4.8 


<-0.6 


< 17.2 


<-1.0 


<28.4 



it seems likely that some absorbers will consist of several com- 
ponents whereas others may be dominated by a single strong 
component. 

The simulations also depend on the scatter in (O V/H I) at 
fixed N(H I), which depends on the scatter in both the metal 
abundances and physical conditions of the absorbers. Using 
C IV absorption, Rauch et al. (1997) and Hellsten et al. (1997) 
find a scatter of around an order of magnitude, while the anal- 
ysis of Dave et al. (1998) suggests that the true scatter may be 
smaller as some of the observed scatter could be attributable to 
fitting uncertainties. The individual systems measured by El- 
lison et al. (2000) indicate an rms scatter of around 0.5 dex. 
We generate simulated data with a Gaussian distribution in 
(O V/H I) with an rms scatter, which we denote as ao, of 0.0, 
0.5, and 1.0 dex. The results of our simulations of sample L for 
different values of b and ao are plotted in Figure 10, compared 
with our observed value. For each line in Figure 10 we have 
computed 100 simulated data sets for each point in steps of 0. 1 
in (O V/H I), and the plotted value represents of the mean of 
these simulations. The ^-parameter makes a small difference, 
but (To has a larger effect. If we assume that ao is likely in 
the range 0.5-1.0 and allowing for the uncertainty in b, we get 
agreement with our observations for (O V/H l)« -1 .7 to -0.8. 

Generating simulated datasets for sample I is not as simple. 
Because we do not have accurately measured column densities 
for these lines, we cannot simply calculate the O V column den- 
sities from the H I column densities of the lines in our sample. 
Instead, we use simulated lists of random lines drawn from a 
column density distribution with [3 = -1.41 for N(H I)< 10 143 
cm -2 and /3 = -1.83 for 10 14 ' 3 < N(H I)< 10 16 2 cm -2 , as shown 
in Figure 1 . We assume a ^-parameter distribution as derived 
by Kirkman & Tytler (1997): a Gaussian distribution in b with 
b = 23 km s" 1 and a^ = 14 km s" 1 , and a minimum b given by 
b min = 14 + 4 x log[V(H I)/12.5] km s _1 . The distribution of b 
is important because we have used a cut in equivalent width to 



define our sample, so we must know the equivalent width of the 
lines in the simulated data. In our full sample, we have 69 Lya 
lines with 10 14 2 <JV(HI)< 10 16 2 cm . In Figure 1 1 we plot a 
histogram of the equivalent width distribution of these lines, as 
well as the average of 1000 simulated lists of Lya lines in this 
column density range. The notable difference between the two 
is the excess of large equivalent width lines, EW > 0.8 A, in the 
real data. These are likely metal-line systems, as discussed ear- 
lier, though they could also be blended lines. The upper-limit 
in Lya equivalent width of 0.6 A for sample I is thus chosen to 
help assure that this sample is relatively free of metal-line sys- 
tems. At lower equivalent widths, it seems that the simulated 
data are a good approximation of the real data. 

In generating our simulated O V data for sample I, we replace 
our sample absorbers with randomly generated line lists. Given 
that we now have lines with a realistic distribution of H I col- 
umn densities, we can calculate the O V column densities and 
generate simulated data as we did for the lower column densi- 
ties. Figure 12 shows the results of the simulations. At these 
larger O V equivalent widths, ao makes less of a difference, but 
the results are quite sensitive to the assumed ^-parameter, since 
for a small ^-parameter the stronger lines will begin to saturate. 
An effective ^-parameter of 10 km s -1 is almost certainly un- 
realistically low for this sample, since C IV absorption in this 
regime often shows structure, but it is useful to consider as an 
extreme case. Depending on the ^-parameter, we find (O V/H I) 
should be in the range -1.6 to -0.6 for sample I. This agrees 
with what we derived for sample L, giving us confidence in our 
analysis and suggesting that (O V/H I) does not vary more than 
an order of magnitude over the range of H I column densities 
that we have considered, from 10 13,2 to^l0 15 cm" 2 . 

5.2. Other Ions 

To this point we have only analyzed O V since it is the only 
ion we detect in samples L and I. However, in §4.4 we placed 
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90% confidence limits on the equivalent widths of absorption 
from other ions. By creating model data as we have done for 
O V A630, we can convert these 90% confidence upper limits 
on the equivalent widths of metal lines of other ions into con- 
straints on the column density ratio of the corresponding ions to 
H I. For all of the ions we assume an rms scatter in (Xi/H I) of 
0.5 dex for the simulated data. Table 6 lists our derived values 
and limits on (X;/H I) for samples L and I. We do not calculate 
values for sample H because our lack of knowledge of the H I 
column densities makes the results very uncertain. 
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FIG. 10. — Expected O V equivalent widths from simulated data for Lya 
absorbers with 10 13 2 < N(H I)< 10 14 2 crrT 2 , as a function of the mean ratio 
(O V/H I). The different curves correspond to different assumptions for the 
i-parameter for O V and the rms scatter in (O V/H I), denoted as a . Our 
measured value is shown as the horizontal solid line, and the Ict limits on the 
measured value are shown as the horizontal dotted lines. 
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FIG. 11. — Solid Line: Distribution of rest-frame equivalent widths of the 69 
Lya absorption lines in our sample with H I column densities > 10 14 ' 2 cm -2 . 
Dotted Line: Average distribution of rest-frame equivalent widths from 1000 
randomly-generated lists of 69 Lya lines with column densities > 10 14 ' 2 cm -2 , 
assuming the column density distribution shown in Figure 1 and a fo-parameter 
distribution as derived by Kirkman & Tytler (1997). 



5.3. The Oxygen Abundance 

Converting a value for the (O V/H I) into an oxygen abun- 
dance (O/H) requires knowledge of the relative ionization cor- 
rection for the ions O V and H I. This is generally ascertained 



by means of photoionization models in which the IGM is illu- 
minated by the processed radiation of QSOs and/or stars. We 
use the photoionization code CLOUDY (version 94.00; Fer- 
land 1996). Because we have not resolved individual O V ab- 
sorption components, in using these models to infer (O/H) we 
must necessarily assume that the absorbers are predominantly 
single-phase; i.e., the O V and H I absorption occurs in the 
same gas. The main inputs to the models are the density and 
temperature of the absorbers as a function of H I column den- 
sity, and the shape and normalization of the ionizing contin- 
uum. For the physical conditions in the absorbers we use the 
results of Hellsten et al. (1998), which are based on hydrody- 
namical simulations as described by Katz, Weinberg, & Hern- 
quist (1996). At z = 3, Hellsten et al. (1998) find that the typ- 
ical density of the absorbers in their model are well-described 
by a simple analytic fit, logn H = -14.8 + 0.71og/V(H I) assum- 
ing Qbh 1 = 0.0125. The typical temperature of the absorbers 
is described by logT = 7.17 + 0.651og«H for log«H < _ 3.8 and 
logT = 3.18-0.41ogn H forlogn H > -3.8. The temperature re- 
lation is somewhat dependent on the ionization history of the 
universe (Hui & Gnedin 1997), particularly for underdense gas, 
but for the absorbers and redshifts we are studying this is a 
small effect and our results are fairly insensitive to the tem- 
perature. 
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FIG. 12. — Same as Figure 10 but for absorbers with N(R I)> 10 14 2 cm -2 
and a Lya equivalent width < 0.6 A. In contrast to the weaker lines in Fig- 
ure 10, the scatter cfq makes only a small difference but the result is more 
sensitive to the assumed ^-parameter, since many of the O V features would 
begin to saturate if b were small. 

For our nominal ionizing spectrum we use the standard 
Haardt & Madau (1996, hereafter HM96) spectrum. Specifi- 
cally, we use their result for z = 2, although there is very little 
evolution in this spectrum between z ~ 2 and 3. HM96 assume 
that QSOs are the source of the ionizing radiation for the IGM. 
They assume a QSO spectral index of «euv = —1-5 {f v oc v a ) 
for the continuum shortward of the Lyman limit, consistent with 
the oeuv = —1.57 ±0.17 for radio-quiet QSOs derived by Telfer 
et al. (2002). Some recent hydrodynamical simulations of the 
IGM, however, suggest that a softer ionizing continuum than 
HM96 is necessary to explain the observed ratio of He II to H I 
in the IGM. Zhang et al. (1997) require less He II ionizing ra- 
diation by a factor of ~ 4, while Theuns et al. (1998) require 
a factor of ~ 2 less. Croft et al. (1997) resolve the discrep- 
ancy by increasing fi^/i 2 by an amount consistent with recent 
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determinations of Qbh? from D/H (O'Meara et al. 2001) and 
the cosmic microwave background (de Bernardis et al. 2002). 
Because the ions of interest to our discussion, in particular O V, 
are created by photons with wavelengths shortward of the He II 
break, our results are sensitive to differences in the He II pho- 
toionization rate. We create different ionizing spectra for our 
photoionization models by first creating a simple analytic fit to 
the HM96 spectrum, which we find produces indistinguishable 
results for the ions of interest. We then create softer spectra by 
simply increasing the strength of the break at the He II edge and 
then extending the softer spectrum to where it meets the HM96 
spectrum in the X-ray regime. The original HM96 spectrum 
and our model spectra are shown in Figure 13. For all these 
models we adopt the normalization of the ionizing continuum 
of J v = 10~ 21,3 ergs cm" 2 s" 1 Hz" 1 sr" 1 at 1 Ryd as derived by 
HM96 and in good agreement with Giallongo et al. (1996). 
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FIG. 13. — Various ionizing spectral shapes considered in this work. The 
solid line shows the spectrum of Haardt & Madau (1996) at z = 2. The dotted 
line shows an approximation to this spectrum using just a few points. These 
two spectra yield virtually indistinguishable results for the ions we consider. 
The dashed lines show the same model spectrum, except we have artificially 
increased the break at the He II ionization edge by factors of two and four. 

In Figure 14 we plot the ionized fractions of O IV, O V, O VI, 
and C IV relative to H I for our photoionization models. For 
the HM96 spectrum, (O V/H l)-(0/H) remains relatively con- 
stant for 10 13 <Af(Hl)< 10 14 8 cm" 2 at around 3.6. Given our 
values of (O V/H I) for sample L, this implies (O/H) around 
-5.3 to -4.4, or [O/H] around -2.2 to -1.3 with respect to the 
standard solar oxygen abundance. In the softer ionizing spectra 
models, slightly less oxygen is required to explain our mea- 
surement for sample L but slightly more oxygen for sample I. 
At larger H I column densities, (O V/H I) varies more dramati- 
cally with N(H I) and the interpretation of our measured value 
of (O V/H I) over a range of A^(H I) is less obvious. 

To include the effect of a changing ionization state with 
N(H I), we create simulations using these specific ionization 
models to determine 7V(0 V) for the simulated absorption where 
the free parameter is now (O/H) rather than (O V/H I) . We cal- 
culate model equivalent widths for O IV A788 as well to com- 
pare with our limits since this can provide an interesting con- 
straint on the ionization state. In Figure 15 we plot the equiva- 
lent widths of O V and O IV as a function of (O/H) for sample 
L, assuming b = 50 km s" 1 and gq = 0.5 dex. We do not plot 
the results for different values of the ^-parameter and <jq, but 



the effects are very similar to those seen in Figure 10; i.e., the 
required (O V/H I) is more by around 0.1 for b = 10 km s" 1 , and 
less by around 0.4 for do = 10 dex. As expected, the required 
(O/H) agrees with what we estimated above, with slightly less 
oxygen required for the softer spectra. For none of these ion- 
izing spectra does the required amount of oxygen imply that 
there should be a clearly observable feature for O IV, although 
for the He II break x 4 spectrum the expected equivalent width 
of O IV A788 is around our 90% confidence limit. This compar- 
ison of O V to O IV is virtually independent of our assumptions 
about the ^-parameter and ao since the O V and O IV equivalent 
widths essentially scale by the same factor for changes in these 
parameters. 
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FIG. 14. — Ionized fractions of various ions relative to H I. The solid lines 
are for O IV-O VI and are labeled in the diagram. The dotted line is for C IV. 
The models correspond to the ionizing spectra shown in Figure 13. A normal- 
ization of J v = 10~ 21 3 ergs cnT 2 s" 1 
models. 



Hz 1 sr 1 at 1 Ryd is used for all three 



The O IV limit is more interesting for sample I. The predicted 
equivalent widths are plotted in Figure 16 assuming <jq = 0.5 
dex and b = 50 km s" 1 , which as we argued before is proba- 
bly more realistic than b = 10 km s" 1 for 7V(H I)> 10 14 ' 2 cm" 2 . 
Again, the effects of changing the ^-parameter and <to are sim- 
ilar to those in Figure 12. For these absorbers, the simulated 
data suggest that the observed lack of O IV favors a hard ion- 
izing spectrum. For the HM96 spectrum the predicted equiva- 
lent width is only slightly higher than our 90% confidence limit. 
However, for the He II break x 4 spectrum, the simulations sug- 
gest an equivalent width of O IV of more than three times our 
limit. 

This constraint is sensitive to the normalization of the ioniz- 
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ing continuum. An increase in the ionizing flux will decrease 
(O IV/O V) resulting in less predicted O IV. We rerun the sim- 
ulations using J v = 10" 21 ergs cm" 2 s" 1 Hz" 1 sr" 1 at 1 Ryd as 
found by Cooke, Espey, & Carswell (1997). The resulting 
(O IV/O V) is smaller by around 30-40%, consistent with the 
HM96 spectrum but still lying above our upper limit for the 
softer spectra. Density gradients in the absorbers can also influ- 
ence the results. By using a constant density, we are assuming 
that the absorption takes place primarily in the densest portion 
of the absorber. This should be a good approximation for a 
given ion provided that the volume density of the ion increases 
with total density. However, sample I contains absorbers with 
H I densities near the point where O V peaks. If a significant 
amount of the absorption were to occur in the density wings of 
the absorber, the observed (O IV/O V) would be less than what 
we have calculated for a constant-density model. 
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FIG. 15.— Expected equivalent widths of O V A630 and O IV A788 as a 
function of (O/H) forLya absorbers with 10 13 2 < N(H I)< 10 14 2 cnT 2 . The 
different lines correspond to different assumed ionization fractions, as shown 
in Figure 14, based on different ionizing spectra. For all the simulated data we 
assume an rms scatter in (O /H) of 0.5 dex and a ^-parameter of 50 km s" 1 . Our 
90% confidence upper limit on the equivalent width of O IV A788 is shown as 
the dotted line in the top panel. Our measured value for O V is shown in the 
bottom panel as it was in Figure 10. 



60 



S 40 

H 

* 20 





100 
80 



6 



> 
o 



60 
40 
20 




— i — 1 — 1 — 1 — i — 1 — 1 — 1 — r 

HM96 

Hell break x 2 

Hell break x 4 




-5.4 



-5.2 



-4.8 

(O/H) 



4.6 



-4.4 



A. 2 



Fig. 16. 
a Lya equivalent width < 0.6 A. 



An additional consideration is our ignorance of the column 
densities of the lines in sample I. Because there are only 49 
lines in the sample, it is possible that the true distribution of 
H I column densities is significantly different than our simu- 
lated lists. In particular, there could be a relative lack in the 
true sample of high-column-density systems where (O IV/O V) 
is highest. Our simulations suggest that the uncertainty in the 
predicted O IV equivalent width due to statistical fluctuations 
in the column-density distribution is comparable to the random 
noise. The additional uncertainly for O V is much smaller since 
in our photoionization models N(0 V) varies much less with 
7V(H I) over the range of H I column densities spanned by sam- 
ple I. Models with predicted O IV equivalent widths only mod- 
erately above our stated limit are therefore consistent with our 
results. However, it is quite difficult to reconcile our limit on 
O IV with the He II break x 4 spectrum. We conclude that our 
data are probably inconsistent with such a soft ionizing spec- 
trum and that the true ionizing background is not significantly 
softer than that of HM96. 

As we noted above, early observational constraints on 
7V(He II) /N(H I), often called 77, led us to consider ionizing con- 
tinua softer than HM96. The strongest constraint on r\ generally 
considered in the literature is from the He II opacity measure- 
ment of Davidsen, Kriss, & Zheng (1996). Recently, Kriss et al. 
(2001) analyzed a high-resolution spectrum of the He II absorb- 
ing region with the Far Ultraviolet Spectroscopic Explorer en- 
abling for the first time the explicit measurement of 77 for indi- 
vidual absorbers. Using both measured values of r\ and lower 
limits for He II absorbers with no corresponding H I absorption, 
they find a mean r\ of <~ 80. However, the mean rj for absorbers 
with detectable H I (7V(H I)> 10 123 cm" 2 ) is ~ 30. Because our 
sample is selected by H I absorption, it is appropriate to com- 
pare to this value. Comparing to the models of Fardal, Giroux, 
& Shull (1998), rj = 30 suggests a QSO source spectrum with 
ckeuv around -1 .6, consistent with Telfer et al. (2002) and simi- 
lar to HM96. Thus the fact that our data support a hard ionizing 
continuum is consistent with the He II results if one takes into 
account that an H I-selected sample preferentially selects ab- 
sorbers with lower values of rj which are likely photoionized by 
hard, QSO-like radiation. 

Including uncertainty in the scatter, our models suggest 
(O/H) of around -5.2 to -4.6 for sample I, although the re- 
quired (O/H) could be much higher if the effective ^-parameter 
is much smaller than 50 km s" 1 . This agrees with the results for 
sample L, for which the results are more robust due to the fact 
that the H I column densities of the constituent absorbers are di- 
rectly measured. We therefore conclude that the typical (O/H) 
for the IGM likely lies in the range -5.3 to -4.4, as derived 
for sample L. The abundance with respect to the standard solar 
value is thus [O/H] around -2.2 to —1.3. 

Having derived this abundance, it is interesting to compare 
this to results from studies of C IV. Using the optical depth ra- 
tio technique, Ellison et al. (2000) find that (C IV/H l)= -2.6 
for T]_ ya ^ 1 provides a good match to their data, in agree- 
ment with what Songaila & Cowie (1996) found for individu- 
ally detectable C IV absorbers. Using our HM96 model for the 
ionization correction for C IV, (C IV/H l}= -2.6 corresponds to 
[C/H] w -2.5. This is in good agreement with what other au- 
thors have derived assuming a HM96 spectrum (Hellsten et al. 
1997; Dave et al. 1998). The inferred [O/C] is in the range 0.3- 
1.2, assuming that the chemical composition does not change 
dramatically from z ~ 2 to 3.5 since our results are for slightly 
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lower redshift than those for C IV. Similar results are obtained 
for the He II break x 2 spectrum. This model implies less car- 
bon, [C/H] « -2.7, to match the C IV results, but also less oxy- 
gen to match our results, resulting in the same [O/C]. Thus, 
since the relative ionization corrections in our models for C IV 
and O V vary little for the H I column densities of interest, 
the uncertainty in [O/C] is dominated by our uncertainty in 
(OV/HI). 

An overabundance of oxygen relative to carbon is expected 
if the IGM is enriched in metals by Type II supernovae. Ed- 
vardsson et al. (1993) find overabundances of oxygen relative 
to carbon by factors of 3-5 in halo stars, consistent with the low 
end of the range we infer. There is also evidence that this over- 
abundance may increase with decreasing metallicity (Israelian, 
Garcia Lopez, & Rebolo 1998). A relative overabundance of 
oxygen is expected from Type II supernova yields (Woosley 
& Weaver 1995). However, large values of [O/C] may re- 
quire nucleosynthesis via pair instability supernovae (Heger & 
Woosley 2002) from very massive stars (m > 100 M@) in the 
early universe (Abia et al. 2001 ; Schaerer 2002). A bimodal ini- 
tial mass function for Populations III stars, including a peak at 
~ 100 Mq, is predicted by the hydrodynamical simulations of 
Nakamura & Umemura (2001). An additional predicted conse- 
quence of nucleosynthesis by very massive stars is a large abun- 
dance of silicon, [Si/C] ~ 1.0. In contrast, Songaila & Cowie 
(1996) infer a typical value of [Si/C] ~ 0.4 for metal-line sys- 
tems from the observed Si IV / C IV ratios. There is evidence 
for absorbers with [Si/C] > 1.0 at high redshift (z > 3), but 
the assumed implausibility of such abundance ratios are used to 
argue instead for a substantial softening of the ionizing back- 
ground (Savaglio et al. 1997; Giroux & Shull 1997; Songaila 
1998). Further exploration of the ionizing background at z > 3, 
such as from measurements of the He II / H I and O IV / O V ra- 
tios to be made possible by the Cosmic Origins Spectrograph, 
will help to clarify this issue. 

Evidence for an overabundance of oxygen at high redshift 
has been found by other authors. Reimers et al. (1992) find 
[O/C] > 0.6 for Lyman-limit systems. Using the optical depth 
technique, Dave et al. (1998) find a good fit for their data us- 
ing [O/C] = 0.5 for C IV absorbers at z > 3. However, based 
on the lack of narrow O VI absorption lines in their data, they 
claim that the overall metallicity must drop dramatically, by a 
factor of at least ~ 3, for absorbers with N(H I)< 10 14 5 cm" 2 . 
In contrast, our results suggest that there is a large abundance of 
metals at A^(H I) < 10 14,2 cm" 2 , consistent with what is found at 
higher H I column densities. The O VI detection of Schaye et al. 
(2000) supports this result by suggesting that the enrichment of 
the IGM does indeed extend to densities below the universal 
mean. 

It is instructive to compare our results to Schaye et al. (2000). 
Assuming (O V/H l}= -1.1 and using the HM96 model for the 
ionization correction, for Lya optical depths ranging from 0.1 
to 1 .0, corresponding to 10 136 < N(U I)< 10 14 6 cm" 2 for b = 30 
kms" 1 , we would predict (OVI/Hl) from around -0.5 for 
T]_ ya = 0.1 to —1.3 for r Lya = 1.0. The corresponding ratio of 
the optical depth of O VI A 1032 absorption to Lya would range 
from —1.1 to —1.9. This agrees well with the apparent optical 
depths of O VI measured by Schaye et al. (2000) for the same 
redshift range, though as they point out, the true optical depth of 
O VI could be considerably different from the apparent optical 
depth due to the effects of noise and Lyman-line contamination. 

The existence of metals at low column densities addresses 



the issue of the source of the metals in the IGM. The enrich- 
ment of the IGM is thought to occur in one of two basic ways: 
(1) by a high-redshift era of small star-forming regions, either 
Population III stars (Ostriker & Gnedin 1996) or protogalaxies 
(Madau, Ferrara, & Rees 2001), or (2) by in-situ star forma- 
tion, either in the Lya absorbers themselves or nearby galaxies. 
The primary distinction between the two scenarios is the vol- 
ume extent of the enrichment. While early supernovae could 
be expected to pre-enrich the universe in a fairly uniform way, 
it is difficult to enrich the most diffuse regions of the IGM by 
in-situ star formation. For example, Gnedin & Ostriker (1997) 
predict a strong dropoff in metallicity at column densities be- 
low around 10 13,5 -10 14,5 cm" 2 at z = 3. Although our sample 
L includes absorbers down to 10 132 cm" 2 , we cannot be sure 
that there are metals at such low column densities, though the 
fact that the S/N of the detection does not decrease until we 
include absorbers with N(H I)< 10 13 2 cm" 2 certainly provides 
evidence that there are metals at low column density. Our re- 
sults apply to slightly lower redshifts than this prediction, so 
there is more time for metals to diffuse into less dense regions, 
and the absorbers that we observe correspond to slightly larger 
overdensities for the same N(H I). Given these considerations, 
we cannot conclusively claim that our results are inconsistent 
with in-situ enrichment scenarios, which are currently an active 
area of study (Ferrara, Pettini, & Shchekinov 2000; Aguirre 
et al. 2001). However, it is fair to say that our detection of 
O V in sample L, in conjuction with previous detections of C IV 
and O VI in similar density regimes, provide evidence for the 
presence of metals at quite low column densities, which would 
favor a uniform enrichment as provided by an early era of star 
formation. 

6. SUMMARY 

We have used HST FOS spectra of four QSOs to search for 
absorption features associated with Lya forest absorbers in the 
EUV in the redshift range 1.6 < z < 2.9. The results of this 
search can be summed up as follows: 

1 . We detect O V A630 over a large range of N(H I). Most 
interestingly, we detect O V in a sample of absorbers 
with 10 13 2 < N(U I) < 10 14 2 cm" 2 with greater than 
99% confidence. 

2. We detect O IV A788, O IV A554, O III A833, and He I 
A584 only for the strongest absorbers in our sample, 
those with Lya equivalent widths > 0.6 A. 

3. We find no evidence for N IV A765, Ne V A568, 

Ne VI A559, Ne VIII AA770, 780, or Mg X AA610, 625 
absorption in any of our samples. 

For absorbers with Lya equivalent widths < 0.6 A, the O V 
detections imply (O V/H l)« —1.7 to -0.6 based on our simu- 
lated data, where this range allows for uncertainties in the as- 
sumed ^-parameter and scatter in (O V/H I). This result implies 
that O V is typically a factor of ~ 10-100 more abundant than 
C IV. The O V A630 line thus has EW/A a factor of - 30-300 
larger than C IV A1548, or ranging in strength from around that 
of LyS to Ly/3, making it an excellent tracer of metal content. 

Using photoionization models to calculate the ionization cor- 
rection, we find that the oxygen abundance in the IGM is 
[O/H] « -2.2 to -1.3, implying [O/C] w 0.3 to 1.2, consistent 
with what is found in halo stars and Lyman-limit systems. The 
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overabundance of oxygen suggests Type II supernova enrich- 
ment, but an unusual stellar initial mass function resulting in a 
significant contribution from pair instability supernovae of very 
massive (m > 100 M Q ) Population III stars may be necessary. 

The fact that we find no evidence for O IV A788 absorption 
except in the strongest systems provides an interesting con- 
straint on the ionizing background spectrum, specifically that 
it is unlikely to be more than a factor of ~ 2 softer than the 
HM96 spectrum at 4 Ryd. A hard ionizing spectrum is consis- 
tent with measurements of He II absorption for absorbers with 
detectable H I. 

We conclude that studying the absorption of the IGM in the 
rest-frame EUV is not only a useful but essential tool for gain- 
ing a more complete understanding of the metal content and 



ionization of the IGM. From longer wavelength data alone one 
cannot place useful constraints on absorption from multiple 
ionization stages of the same element in the diffuse IGM. Ulti- 
mately, a more complete understanding requires such analysis, 
and the clearest path to this end is through the simultaneous 
study of multiple ionization stages of oxygen (O III-O VI). We 
look forward to the installment of the Cosmic Origins Spectro- 
graph on the HST, currently planned for 2004, which will enor- 
mously increase our ability to observe weak oxygen features on 
an individual basis. 



We are grateful to A. Songaila for providing us with the Keck 
spectrum of HE 2347-4342. 
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